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Two-particle interferometry, a second-order interference effect, is explored as another possible tool 
to distinguish between massive Dirac and Majorana neutrinos. A simple theoretical framework is 
discussed in the context of several gedanken experiments. The method can in principle provide both 
the mass scale and the quantum nature of the neutrino for a certain class of incoherent left-handed 
source currents. 
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Two contemporary problems in neutrino physics are 
determining the absolute mass of the neutrino and dis- 
covering if the neutrino is a Dirac or Majorana fermion 
P|. The existence of neutrino mass has been estab- 
lished through oscillation experiments such as Super- 
Kamiokande, SNO, and KamLAND 0, H H, which 
have successfully extracted the differences of the squared 
masses between the energy eigenstates. Various exper- 
imental approaches, such as tritium decay d Q and 
cosmological background studies 0, are capable of ex- 
tracting the kinetic mass of the electron neutrino and 
mass sum of the neutrino energy eigenstates respectively. 
While these experiments have been able to put an ever- 
improving upper limit on the neutrino mass, they provide 
no information about the neutrino's Majorana or Dirac 
nature. One powerful approach currently used to deter- 
mine the quantum nature of the neutrino is neutrinoless 
double beta decay {(3f3{0v)) J2J. The decay rate is pro- 
portional to the effective mass of the neutrino and only 
proceeds if the neutrino is a Majorana particle. Other 
interesting methods to address these fundamental ques- 
tions have been explored since the 1950s but explor- 
ing new ideas may be beneficial. 

This Letter investigates another technique, two- 
particle intensity interferometry, that theoretically pro- 
vides information about the mass and nature of the neu- 
trino. This form of interferometry has been used exten- 
sively in many areas of physics and has served to cross- 
pollinate ideas in different sub-fields for over forty years. 
It is natural to wonder what role this technology might 
play in neutrino physics. 

Intensity interferometry was originally developed by 
Robert Hanbury Brown and Richard Twiss (HBT) as an 
alternative to Michelson interferometry to measure the 
angular sizes of stars in radio astronomy E3- The ideas 
behind intensity interferometry were eventually quantum 
mechanically applied to photons, rather than classical 
radio waves, instigating a revolution in modern quantum 
optics E3 • The technology was independently developed 
in momentum space for final state particles in elementary 
particle physics and is sometimes called femtoscopy in 
that context El El El El El E3 



The basic observation in two-particle interferometry is 
pairs of incoherently generated indistinguishable bosons 
tend to bunch close in phase space while similarly gener- 
ated fermions tend to anti-bunch. What "close" means 
exactly depends on the scale and geometry of the problem 
and in what space one is performing the measurement. 
The sensitivity of the effect to the quantum statistics 
obeyed by measured pairs, in particular the tendency for 
incoherent fermions to anti-bunch in phase space, is of 
interest here in an attempt to determine the quantum 
nature of the neutrino. 

A common physical-observable in intensity interferom- 
etry is the two-particle correlation function, C2, which is 
a measure of the degree of independence between two 
events in some variable of interest, such as momentum, 
space, or time. The two-particle correlation function can 
be written 
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P(l,2) Ti[pa\a\a k a q ] 
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where P(l,2) represents the joint probability of measur- 
ing two events while P(i) represents the individual prob- 
abilities of events i = 1, 2 and can be naturally general- 
ized to higher order correlations. The explicit momen- 
tum space form of C2 on the right hand side of Eq. Q 
highlights the basic physical components of the correla- 
tion function. Tacitly contained in the density matrix, 
p, when projected as a Wigner function, are the space- 
time geometry of the source, the source dynamics, and 
any pairwise interactions. The quantum statistics of the 
particles are determined by the (anti-)commutation rela- 
tions of the creation and annihilation operators, a' and 
a. When normalized to the single particle distributions 
as shown, Ci is proportional to the relative probability 
for a joint two-particle measurement as compared to two 
single-particle measurements. If the measurements are 
independent, then 6*2 = 1 . If the measurements are corre- 
lated, C2 deviates from unity. 

As Eq. implies, there are many possible approaches 
one can use to obtain an explicit expression for the cor- 
relation function. A particularly simple form for Eq. {Q 
that illustrates the essential physics is given by the 
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Koonin-Pratt equation [l^.lT^| 

C 2 = J tfR^x^x^piR). (2) 

The equation assumes an incoherent emission of parti- 
cles from a normalized pair distribution p(R) where R 
is the vector separation between the pairs at the source. 
For simplicity, time has been implicitly integrated out of 
Eq. However, the formalism can be expanded to in- 
clude correlations in the time domain. The two-particle 
wave function, ip(xi, X2), contains information about the 
quantum statistics and any pairwise interactions. Work- 
ing in natural units (c = % = 1), if a pair of free iden- 
tical fermions in any specific triplet spin configuration is 
considered, the spatial part of the wave function will be 
antisymmetric upon label exchange (assuming no other 
quantum numbers are involved) and given by the usual 
plane wave solution 

^(xi,f 2 ) = J^e^^V*-* 1 - e&'^'e**-* 1 ). (3) 

One interprets this two-particle wave function as the am- 
plitude for particles emitted at points x\ and x 2 to be 
measured with momenta p a and pb- For free particles, 
C 2 is simply related to the cosine transform of the inco- 
herent pairwise source distribution, p{R). 

If two identical free fermions are emitted from exactly 
two point sources separated by R, Eq. J3J) can be written 

C 2 (Q) = l-^cos(Q-fe) (4) 

where Q = p a — Pb and 5x = x\ — x 2 . The parameter 
£ = 1 for triplet spin states and —1 for singlet states. 
If the system is spin-averaged, then £ = |. Notice in 
the triplet case C 2 (Q = 0) = and the fermions are 
anti-correlated if in the same momentum state. Because 
the emission is incoherent, and there are no interactions, 
the correlations arise only from the quantum statistics 
obeyed by the particles. The scale of the correlation is 
set by the source size. It is instructive to note that for 
non-identical particles, where the wave function has no 
particular symmetry, C 2 — 1 for all Q. 

Let's examine a useful limit of Eq. (0J that will be 
used later for a series of gedanken experiments. Consider 
two point sources of fermions separated by a distance R 
and measured by a pair of distant detectors separated by 
d. The source and detector are a distance L from each 
other such that L ^ R 3> d. That is, there are well- 
separated sources far away from a relatively close pair 
of detectors. Assume a pair of single-mode fermions. In 
this limit Eq. becomes 

C 2 (d) = l-£cos(A0d/A). (5) 

Although the particles here are fermions, this is similar 
to the original HBT experiment used to measure the an- 
gular size of stars. The correlation function is measured 



at different detector separations, d, for waves of known 
wavelength, A. From the shape of C 2 (d), the angular size, 
A9, can be extracted. 

Imagine not knowing a priori the quantum nature of 
the particles being measured, but instead knowing some 
other information such as the angle subtended by the 
source relative to the detectors. In that case, using 
Eq. 0, one would fix the angular size and wavelength 
but then look at the behavior of C 2 as the distance be- 
tween detectors approached zero to determine the quan- 
tum statistics obeyed by the particles of interest. 

Can two-particle interferometry be applied to neutri- 
nos to determine if they are Dirac of Majorana parti- 
cles? Let's examine four variations of a simple gedanken 
experiment, labeled A through D below, to answer this 
question. A summary of the relevant formulae and the 
ability of the four cases to resolve the neutrino mass and 
nature are outlined in Table [I] For simplicity, only one 
neutrino flavor with one mass eigenstate is considered 
and oscillations are ignored. 

It will be helpful to remember for the cases below that 
while Majorana neutrinos are their own antiparticle (the 
field operators transform to themselves under a charge 
conjugation operation), the left-handed weak source cur- 
rents creating them will generate final state particles with 
a handedness as if they were Dirac fermions 0. 

TABLE I: The two-particle correlation function for Dirac, 
C^"(d), and Majorana, C 2 a ^(d), neutrinos are shown for var- 
ious situations. Where £ alone is quoted, use Eq. @. An 
entry of C2 = 1 indicates no correlation. The helicity column 
notes if detectors are filtering on same, opposite, or averaged 
final state helicities. The final rows provide an overview of 
the case-by-case physics capability to determine the neutrino 
mass or discover the neutrino nature. Case A: m — 0, iden- 
tical sources; Case B: m = 0, distinguishable sources; Case 
C: m / 0, identical sources; Case D: m 7^ 0, distinguishable 
sources. See the text for a detailed case-by-case discussion. 
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First, in case A, consider a massless neutrino and a ge- 
ometric setup like that describing Eq. © : well-separated 
sources far away from close detectors. Imagine two reac- 
tors acting as incoherent point sources of indistinguish- 
able particles normally called Dirac antineutrinos. Two 



3 



distant detectors are separated by a distance d. Relative 
to the detectors, the reactor pair subtends a known angle 
AO. There are two situations: one where the sources emit 
right-handed massless Dirac antincutrinos and another 
where the sources emit right-handed massless Majorana 
neutrinos. In this case, measuring C2 cannot distinguish 
between Dirac and Majorana particles. The measured 
correlation function will be equal to that in Eq. Q with 
£ = 1 and will give the same result for both the Dirac 
and Majorana cases. This is because quantum indistin- 
guishability applies equally well for the two situations 
and the two-particle wave function will be identical in 
both cases. Indeed, this is a sanity check because in 
the massless limit, Dirac and Majorana particles cannot 
be distinguished based on the Practical Dirac-Majorana 
Confusion Theorem |igj . 

Next, for case B, consider massless neutrinos with a 
similar geometric source-detector setup as above except 
with one of the reactor sources being replaced by a "small 
sun" . That is, there are two sources emitting distinguish- 
able objects: one, an incoherent point source of particles 
normally called Dirac neutrinos, and another that would 
again be Dirac antineutrinos. But similar to case A, not 
knowing the neutrino nature, there is no way to use C2 
to determine if there is one source of Dirac neutrinos 
and another of Dirac antineutrinos or if there is a pair of 
sources emitting Majorana neutrinos of opposite hand- 
edness. The correlation function 02(d) = 1 for both sce- 
narios. This is because the two-particle wave function 
for either has no special symmetry. That is, it factorizes 
and the particles are not entangled at the detector. From 
Eq. (0), if the normalized wave function factorizes, the 
correlation function becomes unity. 

For C and D let's consider the above two cases again 
but this time give the neutrino a mass that is small com- 
pared to its energy. The presence of mass complicates the 
situation because chirality ("handedness"), is no longer 
the same as helicity. Also, for a realistic Majorana mass 
term, like that introduced in the seesaw mechanism, the 
mass-degenerate four-component Dirac spinor splits into 
two two-component Majorana spinors. For the Majorana 
cases below, we can imagine taking the light doublet, 
keeping in mind that the value of the mass, to, will be 
different in the Dirac and Majorana cases but both will 
still be light compared to the mass of other leptons. 

The primary effect of interest is that left-handed weak 
source currents can now create massive neutrinos and 
antineutrinos of the "wrong" helicity with an amplitude 
that goes like m/E when m <C E. This will determine 
the the probability of measuring an indistinguishable pair 
in the final state that will be treated separately for the 
Dirac and Majorana sources. If the measured fermion 
pair is indistinguishable, the wavefunction must be anti- 
symmetric. The probability of this to occur, which will 
be related to the mass, will determine the strength of 
the two-particle correlation function. This is similar in 



spirit to the considerations in neutral kaon femtoscopy, 
although using different sources, quantum numbers, and 
statistics [2£j. 

With this in mind, consider case C where the source- 
detector geometry with two reactors is the same as case 
A. However, this time each reactor is the source of ei- 
ther Dirac antineutrinos of mixed helicity or Majorana 
particles of mixed helicity. For m< £, the helicity mix- 
ture will be mostly A = +1 with some A = —1 in both 
the Dirac and Majorana cases. For this exercise, con- 
sider ideal detectors that are capable of filtering on the 
neutrino helicity. If the detectors filter on identical helic- 
ities in the final state, C2 will be Eq. with £ = 1, the 
same as case A. Particles of opposite helicity are quantum 
mechanically distinguishable, so if the detectors filter on 
opposite helicities then C2 = 1, as in case B. However, if 
the detectors helicity-average particles in the final state, 
the mixed helicty of the source has the effect of intro- 
ducting a helicity "contamination" at the detector and 
there will be quantum distinguishably for a small frac- 
tion of the measurements. This contamination will have 
the effect of diluting the correlation function by a fac- 
tor d(m 2 /E 2 ); so use Eq. © but with f ~ (1 - m 2 /E 2 ) 
for both Dirac and Majorana particles. Again, Dirac 
and Majorana neutrinos cannot be distinguished, but a 
careful helicity-averaged measurement of C2 (d) could, in 
principle, extract the mass by measuring the strength of 
this weak anti-correlation. 

Finally, in case D, revisit the non- identical sources of 
"sun- reactor" geometry of case B, but extend it to the 
massive neutrino case. Because of helicity mixing, the 
quantum distinguishability arguments are similar to C 
but now there are more combinatorics for the Dirac par- 
ticles because of the extra lepton quantum number. Nev- 
ertheless, like the massless case, the Dirac particles are 
always distinguishable at the detector either by helic- 
ity or by lepton number. No matter how one filters on 
the final state, the Dirac particles are distinguishable so 
C 2 = 1. 

If the neutrino is a Majorana particle, however, case 
D will be different. The reactor source will be emitting 
primarily Majorana neutrinos with A = +1 with a small 
component of A = — 1. The sun source will be emitting 
Majorana neutrinos of the opposite degree of contamina- 
tion: mostly A = — 1 with a small A = +1 mixture. Here, 
because the Majorana neutrino is its own antiparticle, all 
emitted neutrinos are just various helicity states of the 
same particle. With a judicious choice of filtering at the 
detector, one could detect a distinct signal compared to 
the Dirac case. For example, if the detectors filter on 
opposite final state helicity, 02(d) = 1 because the parti- 
cles are distinguishable. But if the detectors filter on the 
same helicity, C2 becomes Eq. (JjJJ with £ = 1 . If a helicty- 
averaging is performed in the final state, this introduces 
contamination (more severe than case C) that will reduce 
the correlation strength. The probability of measuring 
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two equal helicity states with open final-state helicity fil- 
ters scales like m 2 /E 2 so use Eq. J5J with £ ~ m 2 /E 2 . 
The neutrinos would very nearly be anti-correlated at 
small d, with only a slight deviation given by £. 

Let's entertain some experimental considerations. The 
primary concerns are data rate, detector efficiency, and 
energy resolution. The above discussion assumed infi- 
nite energy resolution to resolve neutrinos of an arbitrary 
wavelength with no loss of fidelity or smearing. This as- 
sumption, using Heisenberg's Uncertainty Principle, per- 
mits infinitely slow counting statistics, allowing quantum 
mechanically coherent data to arrive over infinitely long 
time scales. This is clearly an unrealistic practical as- 
sumption. 

The data rates for current experiments such as Kam- 
LAND and SNO are about one event per day. To perform 
the measurement, even assuming copious statistics, the 
ability to measure neutrinos of arbitrary energy, and very 
fine vertex resolution, experiments would require an un- 
physical energy resolution to see the effect as described. 
Conversely, using AEAt ~ (eV)(fs) it can be seen that 
even with extremely good, but still physical, energy reso- 
lutions, eV or keV, an experiment needs to measure neu- 
trino pairs separated by times on the order of femto- 
to attoseconds - a rate approaching weak-charge Am- 
peres of neutrinos. If neutrinos could, in principle, be 
measured experimentally with such copiousness and effi- 
ciency, other methods would mostly likely provide a more 
straightforward path to revealing the neutrino's currently 
unknown properties to the same order in m 2 / E 2 . 

The femtoscopic limit of Eq. J3J [L ^> d 3> R) can also 
be considered. In that limit, neutrinos and antineutrinos 
could be generated from very small sources like those cre- 
ated in a high energy physics collisions. In order to image 
femtometer-sized sources, an experiment would construct 
Ci in momentum space, measuring two or more identified 
or reconstructed inclusive neutrinos per event with a mo- 
mentum resolution of roughly MeV. Finally, the method 
could be applied as an anti-bunching counting experi- 
ment in the time domain, similar to what is done in 
quantum optics with photons. This could be performed 
on a beam of neutrinos and/or antineutrinos, mirroring 
cases A-D above. High flux neutrino-antineutrino beams, 
like those expected from muon colliders, and exceptional 
detection time resolution would be required. 

Based on the gedanken experiments, in particular case 
D, and reviewing Table [i] there is the rather promis- 
ing theoretical result that, with the correct sources and 
filters, two-particle interferometry can obtain both the 



mass and the nature of the neutrino of any flavor us- 
ing a single physical-observable, C%. While the above 
experimental discussion is not meant to be exhaustive, 
it appears the practical requirements currently render 
the method prohibitive and would require a fundamental 
shift in the way neutrinos are detected. 

The author would like to especially thank N. Xu for in- 
spiring and invaluable input. Thanks also to J. L. Klay 
and my other colleagues at LBNL and LLNL for stimulat- 
ing discussions. This work was supported by the Direc- 
tor, Office of Science, of the U.S. Department of Energy 
under Contract No. DE-AC02-05CH11231. 



[1 

[2 

[3 

[4 

[5 
[6 

V. 

[§ 

[9 
[10 

[11 
[12 

[13 
[14 
[15, 
[16 

[17; 
[is; 

[19 
[20 



* Electronic address: tdgutierrez@lbl.gov 
] S. J. Freedman and B. Kayser, Joint Study on the Fu- 
ture of Neutrino Physics: The Neutrino Matrix, (APS 
Multidivisional Neutrino Study, 2005). 
] Y. Fukuda et al. (Super-Kamiokande), Phys. Rev. Lett. 

81, 1562 (1998), hep-ex/9807003. 
.] Q. R. Ahmad et al. (SNO), Phys. Rev. Lett. 87, 071301 

(2001), nucl-ex/0106015. 
] K. Eguchi et al. (KamLAND), Phys. Rev. Lett. 90, 

021802 (2003), hep-ex/0212021. 
.] V. M. Lobashev et al., Phys. Lett. B460, 227 (1999). 
■] C. Kraus et al., Eur. Phys. J. C40, 447 (2005), hep- 
ex/0412056. 

'] S. Hannestad, Phys. Rev. D66, 125011 (2002), astro- 
ph/0205223. 

■] B. Kayser, The Physics of Massive Neutrinos (World Sci- 
entific, 1989). 
'] R. Davis, Phys. Rev. 97, 766 (1955). 
i] R. Hanbury Brown and R. Q. Twiss, Nature 178, 1046 
(1956). 

] R. J. Glauber, Phys. Rev. 130, 2529 (1963). 
] G. Goldhaber, S. Goldhaber, W.-Y. Lee, and A. Pais, 

Phys. Rev. 120, 300 (1960). 
.] D. H. Boal, C. K. Gelbke, and B. K. Jennings, Rev. Mod. 

Phys. 62, 553 (1990). 
:] G. Baym, Acta Phys. Polon. B29, 1839 (1998), nucl- 
th/9804026. 

] M. Gyulassy, S. K. KaufTmann, and L. W. Wilson, Phys. 

Rev. C20, 2267 (1979). 
■] U. W. Heinz and B. V. Jacak, Ann. Rev. Nucl. Part. Sci. 

49, 529 (1999). 
'] M. A. Lisa, S. Pratt, R. Soltz, and U. Wiedemann, Ann. 

Rev. Nucl. Part. Sci. 55, 357 (2005), nucl-ex/0505014. 
.] D. A. Brown and P. Danielewicz, Phys. Lett. B398, 252 

(1997), nucl-th/9701010. 
»] B. Kayser, Phys. Rev. D26, 1662 (1982). 
'] S. Bekele, J. Phys. G30, S229 (2004). 



